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ABSTRACT 

The  contact  problem  associated  with  a bolt  bearing  on  a circular  fric- 
tionless hole  in  a rectangular  orthotropic  plate  has  been  examined.  The 
problem  has  been  formulated  by  the  well-known  methods  of  Muskhelishvili  and 
Lekhnitskii  using  functions  of  complex  variables.  The  orthotropic  stress 
functions  are  expanded  in  Laurent  series  and  boundary  conditions  satisfied 
by  boundary  collocation.  A simple  iterative  process  has  been  utilized  with 
success  to  find  the  contact  angle.  The  process  utilizes  the  physical  fact 
that  the  radial  stress  in  the  contact  region  is  compressive  and  reduces  to 
zero  in  a well-behaved  manner  at  the  point  where  the  bolt  begins  to  sepa- 
rate from  the  plate.  The  problem  was  motivated  by  the  growing  need  in 
composite  designs  to  gain  analytical  understanding  of  the  joint.  Some 
interesting  results  concerning  the  stress  concentrations  occurring  in 
plates  with  various  geometries  and  material  parameters  are  presented.  The 
results  show  the  importance  of  complex  interactions  and  the  necessity  of 
employing  the  particular  method  to  estimate  stress  concentrations  in  prob- 
lems of  this  nature. 
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Poisson’s  ratio 

INTRODUCTION 


With  the  advent  of  modern  composites  it  has  become  important  to  have  analyt- 
ical solutions  which  can  describe  the  elastic  behavior  of  practical  types  of 
composite  structures  under  load.  As  composites  are  orthotropic  in  behavior,  an 
anisotropic  formulation  becomes  imperative.  Earlier  works  done  for  isotropic 
materials  for  a single  bolt  in  an  infinite  plate  by  Bickley1  and  extensions  of 
his  analysis  to  treat  finite  boundaries.2”4  The  present  paper  treats  a mixed- 
boundary value  problem  connected  with  the  rigid  bolt  in  a finite  rectangular  plate 
using  the  Lekhnitskii  complex  variable  approach.5  For  plane  problems  of  orthotropy, 
Lekhnitskii5  has  developed  a solution  in  terms  of  analytic  functions  of  complex 
variable  theory  along  the  lines  of  the  well-known  techniques  of  Muskhelishvili . 5 

The  problem  considered  here  is  the  problem  of  a rigid  bolt  in  a finite  rect- 
angular sheet.  The  load  is  imposed  by  the  bolt  and  its  effect  is  represented  as 
a known  displacement  on  the  hole  of  the  same  diameter  as  the  bolt.  The  plate  is 
kept  in  equilibrium  by  a uniform  stress  on  one  edge  of  the  plate.  Two  types  of 
boundary  conditions  have  been  considered  at  the  lateral  outer  boundaries,  one 
corresponding  to  a single  bolt  in  a rectangular  sheet  and  the  other  to  a series 
of  bolts  equally  spaced  in  a direction  normal  to  the  direction  of  the  applied 
load.  The  boundary  conditions  are  satisfied  by  a collocation  method  used  by  the 
author  in  Reference  7. 


FORMULATION 

Standard  elements  of  a mechanical  fastening  system  are  shown  in  Figure  1. 
The  essential  part  is  an  orthotropic  plate  with  a circular  cutout  to  accommodate 
a rigid  circular  bolt  of  the  same  radius.  Force  is  applied  to  the  plate  by  a 
translation  displacement  applied  to  the  bolt,  which  causes  a deformation  on  the 
circular  hole  whose  pattern  is  known.  A frictionless  contact  has  been  assumed. 
The  plate  is  kept  in  equilibrium  by  a uniform  tensile  stress  a ^ applied  on  the 
left  edge  of  the  plate  in  Figure  1.  In  practice  such  plates  are  fabricated  of 
several  layers,  each  having  orthotropic  properties  and  varying  angles  of  incli- 
nation. The  axis  of  symmetry  in  this  particular  problem  is  the  x-axis. 

In  Figure  1 two  types  of  joint  configurations  have  been  shown,  a single  pin 
and  a periodic  array  of  pins.  Figure  lc  shows  an  isolated  pin  well  removed  from 
the  boundaries  where  an  infinite  plate  solution  will  apply.* 


*OPLINGER,  D.  W.,  and  GANDHI,  K.  R.  Analytical  Failure  Studies  in  Mechanically  Fastened.  Composite  Plates.  Unpublished  research. 
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* 


(C)  Pin  in  Infinite  Plate 


Figure  1.  Mechanical  joint  configurations  under  consideration. 


The  problem  has  been  treated  as  a plane  problem  ignoring  the  effect  of  out- 
of-plane  and  shear  deformations.  Mathematically,  the  problem  has  been  formulated 
using  the  functions  of  complex  variables.  Adopting  the  conventional  notation, 
the  material  properties  can  be  described  by  the  elastic  constants  Ej,  E2,  v12, 
v2i,  and  G 1 2 • By  symmetry  of  the  constitutive  matrix,  v12/Ei  = ^21/^2- 

The  formulation  of  the  analysis  will  now  be  briefly  summarized.  The  Airy 
stress  function  for  two-dimensional  problems  of  rectilinear  anisotropy  can  be 
written  vas 

£ 

U(x,y)  = Re  Ux { (zx)  + U2(z2)},  (1) 

where  Ui  and  U2  are  analytic  functions  of  the  complex  variables  and  z2.  The 
variables  z 1 and  z2  correspond  to 

z i = x + siy 

z2  = X + s2y,  (2) 

where  si  and  s2  are  the  complex  roots  of  the  characteristic  equation 

an)j4  - 2a16y3  + (2a12  + a66)  y2  - 2a26y  + a22  = 0.  (3) 

Equation  3 defines  the  condition  required  for  the  Airy  function  U to  satisfy  the 
two-dimensional  compatibility  equation . 

In  Equation  3 the  a^j  are  the  compliance  coefficients  occurring  in  the  generalized 
Hooke ?s  Law.  Assuming  unequal  roots  in  Equation  3,  we  have 

S1  = Ul  = a{  + i 3 i 
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s2  = y2  = “2  + i&2 

s3  = §i  s4  = s2. 

(4) 

I ! 

where  cu  and  3j  are  real  constants.  Without  loss  of  generality  we  assume  that 
$1  > 0 and  32  > 0. 

To  simplify  notation  the  new  functions  introduced  are: 

$1  (zl)  = dUi/dz^,  <j>2  (^2)  = dU^/ dZ2 • (5) 

Then,  from  the  definition  of  the  Airy  stress  function,  the  stresses  can  be 
written  as 

^ | ^ | 

ax  = 2Re  {si<f>i(Zi)  + S2<j>2(z2)} 
oy  = 2 Re  { <f>  1 Cz x !)  + 4>2(z2)^ 

TXy  = -2Re  {si<f>{  (Zj)  + S2<f>2(z2)^>  (6) 

where,  henceforth,  primes  will  be  used  to  denote  differentiation  with  respect  to 
the  indicated  arguments . 

Another  useful  relation  is  the  following  expression  for  the  force  resultant 
as  a function  of  arc  length  on  an  arbitrary  curve: 

(1  + isi) <j>i  (zj)  + (l  + is2) <j>2 (z2)  + (l  + isi)<j>i  (zi)  + (l  + i§2) <j>2  (z2) 

= i/s(Xn+iYn)ds  = fi(s)  + if2(s).  (7) 

The  displacements  are  expressed  as 

u = 2Re[pi<j>1+p2<j>2] 
v = 2Re[qi4>i+q2<j>2]  , 
with  p^  and  q^  (k=l,2)  given  by 

Pk  = ailSk2  + al 2-al 6sk 

qk  = allsk2  + a22/ sk-a26 • (8) 


METHOD  OF  ANALYSIS 

We  now  need  to  determine  <|>i  Cz  i)  and  <f)2(z2)  so  as  to  satisfy  the  stress  and 
displacement  boundary  conditions.  Since  the  circular  boundary  is  transformed  to 
an  elliptical  shape  in  the  and  z2  planes,  the  first  step  consists  of  finding  a 
simple  mapping  function  which  maps  elliptical  contours  in  the  plane  on  a unit 
circle  in  the  plane  of  an  auxiliary  parameter  . The  following  accomplishes 

this  step: 

Zk  = “^k3  = (R(1-iak)/2)Ck  + (R/2)(l+isk)(l/ck)  (9) 
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where  R is  the  radius  of  the  hole.  The  following  are  useful  notations: 


<h(zl)  = <f>ll>CSi)]  •-  <hUl) 
4’2(z2)  = 4>2  C?2)  ] = <f)2  (£2)  • 


(10) 


Thus 


4* i Cz  1 D = 4>i  CCiD/^1'  CCiD 
$2  Cz 2)  = ^2 (?2) /u ' C?2^ • 


(ID 


The  transformed  expression  for  the  stresses  can  now  be  written  as 


0^  = 2Re{  Cs i4> j (Ci)/a)'  (Ci))  + (s2<j>2  (42)/“'  (£2) ) ) 


(12) 


The  resultant  force.  Equation  7,  becomes 


(l  + isi)<|>i  (?i)  + (l+is2)  cf>2  (C2)  + (l+lsi)  cf>i  (?i) 


+ (l+is2) 42 (?2)  ~ ^l(s)  + lf2(s)* 


(13) 


The  boundary  conditions  involved  in  the  analysis  included  specifications  of 
normal  and  tangential  stresses  at  the  rectangular  outer  boundary  together  with 
mixed  boundary  conditions  on  the  circular  inner  boundary.  On  the  inner  boundary 
the  conditions  are: 


0<9<q;u=6  cos0. 

The  second  and  third  of  these  conditions  imply  that  the  contact  angle  is  identified 
with  n.  With  6 taken  to  be  the  rigid  body  displacement  of  the  pin,  the  radial 
displacement  uR  is  given  to  a good  approximation  (neglecting  higher  order  terms 
in  <5/R)  by  6 cos0  in  the  contact  region.  The  value  of  n is  an  unknown  in  the 
problem  and  was  determined  by  a method  described  later. 

On  the  outer  boundary  the  boundary  conditions  are: 


r = R,  0 < 0 < it;  t = 0 (frictionless  case) 

R0 


n<0<i;;aR=O 


x=e,  0 < y < s;  ax  = xxy  = 0 


x = -L,  0 < y < s;  ax  = aA,  xxy  = 0 

-L  < x < e,  y = 0;  v = xxy  = 0 

-L  < x < e,  y = s;  xxy  = 0,  (single  pin)  ay  = 0 

-L  < x < e,  y=s;  xxy  = 0,  (multipin)  = constant. 
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This  implies  that  the  loaded  end  of  the  plate  is  located  at  x = -L.  The  third 
equation  represents  the  assumption  that  the  x-axis  is  a line  of  symmetry  for  the 
problem.  In  the  last  equation  the  single  pin  case  corresponds  to  a free  boundary 
at  y = s,  while  in  the  multipin  (i.e.,  periodic  array)  case  y = s represents  a 
second  line  of  symmetry.  In  the  case  of  the  contact  angle,  the  value  of  v in  the 
latter  situation  is  an  unknown.  More  precisely,  it  is  required  that  the  constant 
value  of  v at  y = s be  adjusted  so  that  the  mean  value  of  Oy  along  this  boundary  is 
zero.  In  the  approach  used  here,  two  trial  values  for  v were  assumed  and  the  total 
vertical  force  fy  on  the  boundary  was  calculated  as  part  of  the  solution.  Since  f 
is  a linear  function  of  v,  a straight  line  relation  was  thereby  established  which 
determined  the  value  of  v required  to  set  fy  (and  therefore  the  mean  Cy)  equal  to 
zero.  In  the  infinite  plate  case,  of  course,  the  outer  boundary  conditions  do  not 
enter.  Correspondingly,  the  positive  powers  in  the  Laurent  series,  Equation  14, 
are  eliminated.  Other  simplifications  which  were  used  to  advantage  in  the  infinite 
plate  case  are  described  in  Equation  13. 

The  boundary  condition  equations  contain  the  unknown  stress  functions  4> l Cc  1 ) 
and  $2  it 2)  > which  are  the  analytic  functions  of  their  arguments  in  a doubly  con- 
nected region.  Since  we  are  dealing  with  this  type  of  region,  a Laurent  series  of 
the  following  form  has  been  selected: 

+00  n 

4>iUl)  = Ai  log  Cl  + ^__0O  (an+i$n)Ci 

+OO  4 ft 

^2(^2)  = a2  1o5  £2  + £=_oo  ^n+^n^2* 

A further  simplification  can  be  made  knowing  the  symmetries  in  the  elastic  proper- 
ties of  composite  materials.  The  following  major  kinds  are  considered  here.  Both 
of  them  ensure  that  x-axis  is  the  axis  of  symmetry: 

a.  when  the  characteristic  roots  sx  and  S2  are  purely  imaginary,  i.e., 
oq  3 CX2  = 0;  and 

b.  when  the  roots  are  related  by  the  equation  S2  “ -§1. 

The  above  symmetries  are  obtained  depending  upon  the  manner  in  which  composite 
materials  are  put  together.  Most  of  the  composite  designs  would  fall  into  one  or 
the  other  category.  In  the  symmetry  of  type  a it  can  be  shown  that  under  that 
special  condition,  an  = yn  = 0,  while  in  the  symmetry  of  type  b,  6n  = -yn  and  Bn  = 
6n.  Further,  the  complex  coefficients  of  the  log  terms  Ax  and  A2  are  determined  by 
imposing  the  single  value  condition  on  fx  and  f2  and  the  displacements  u and  v. 

As  derived  by  Lekhnitskii , 4 the  following  equations  are  relevant.  The  values  of 
Ax  and  A2  depend  entirely  on  the  magnitude  of  the  pin  load  P = 2sa^. 

Ax  + A 2 - Ax  - A 2 = 0 

sxAx  + S2A2  - s xAx  - S2A2  = -P/2TTi 

P 1 Ai  + P 2A2  - PxAl  ” P 2A2  “ 0 

_ & 

qiAx  + q2^2  - qi^x  - q2^2  = 0 

Now  the  problem  simplifies  to  selecting  the  unknowns  an,  3n,  yn,  and  6n  in  Equation 
14  so  that  the  external  boundary  conditions  are  satisfied  with  sufficient  accuracy 
A least-squares  collocation  procedure7  was  used  and  found  satisfactory.  It  consists 
of  the  following  steps. 

f 

1 


p 
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1.  Truncation  of  the  infinite  expansion  in  Equation  14  to  a suitable  finite 
number  of  positive  and  negative  powers. 

2.  Selection  of  points  or  stations  around  the  external  boundary  in  a fairly 
uniform  manner  and  writing  force  equations  at  each  location.  The  number  of  such 
equations  was  always  arranged  to  provide  more  than  twice  the  number  of  unknowns. 

3.  Solutions  of  unknowns  in  such  a manner  that  the  equations  in  step  2 are 
satisfied  in  a least-squares  sense. 


ANGLE  OF  CONTACT 

As  in  any  contact  problem,  the  determination  of  the  region  of  contact  is  a 
part  of  the  solution.  As  the  contact  length  depends  on  the  load,  the  problem 
strictly  speaking  is  nonlinear.  In  our  problem  we  employed  a simple  iterative 
procedure  which,  within  a few  trials,  succeeded  in  arriving  at  a close  approxima- 
tion to  the  true  contact  angle  n starting  with  an  assumed  value  close  to  90°.  The 
method  relies  on  the  obvious  physical  fact  that  the  radial  stress  oT  in  the  contact 
region  is  always  compressive  and  approaches  zero  at  the  end  of  contact  in  a well- 
behaved  manner.  Therefore,  the  procedure  consists  of  starting  with  a value  of  n 
denoted  by  m*  solving  the  problem,  and  examining  the  sign  of  oT  at  0 = rq.  If 
the  sign  of  ar  is  compressive,  it  can  be  concluded  that  the  true  value  of  n is 
greater  than  the  assumed  value  m«  In  either  case  a solution  is  obtained  with  a 
new  trial  value  of  n,  say  ri2*  In  aH  the  cases  tried,  the  two  trial  values,  m 
and  x)2>  were  sufficient  to  give  a relatively  accurate  estimate  of  the  true  value 
of  n-  This  approach  amounts  to  a procedure  similar  to  the  Newton  Raphson  method 
for  nonlinear  algebraic  equations.  In  the  region  close  to  the  true  value  of  r\, 
the  assumed  contact  angle  plotted  against  ar  at  9 = n gives  a nearly  linear  plot. 
The  iterative  method  converges  quite  rapidly. 


NUMERICAL  RESULTS 

On  the  basis  of  the  foregoing  analysis,  a computer  code  was  developed  to  de- 
termine the  coefficients  of  the  stress  functions.  The  procedure  consisted  of 
selecting  discrete  stations  around  the  boundary  at  which  the  prescribed  boundary 
conditions  were  specified.  At  each  of  these  stations  two  equations  expressing  the 
pertinent  boundary  conditions  were  written.  The  stations  were  roughly  equally 
spaced  and  their  numbers  so  arranged  that  they  gave  a large  number  of  equations, 
more  than  twice  the  number  of  known  coefficients  in  the  stress  function.  The  re- 
sulting redundant  system  was  solved  by  the  least-squares  collaction  procedure.7  To 
investigate  the  accuracy  of  the  solution,  some  arbitrary  points  on  the  boundary 
were  examined  for  boundary  conditions.  For  the  range  of  the  results  presented 
encompassing  different  material  and  geometrical  configurations,  a fairly  high  order 
of  accuracy  was  obtained  even  with  the  truncation  of  the  series  to  -10  < n < + 10. 

The  practical  motivation  in  carrying  out  this  analysis  was  to  develop  a tool 
to  predict  the  complicated  behavior  of  mechanical  joints  of  laminated  composites 
typical  of  aircraft  structural  applications.  The  elastic  properties  of  layers  are 
given  in  Table  1.  The  fabrication  of  laminates  for  a joint  is  usually  done  by 
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Table  1.  ASSUMED  MECHANICAL  CONSTANTS  OF 
FIBER-REINFORCED  MATERIALS 


Elastic 

Constants 

Graphite 

Epoxy 

Glass 

Epoxy 

EL’ 

106  psi  (GN/m2) 

25.0  (172) 

7.00  (48.2) 

Er 

106  psi  (GN/m2) 

1.0  (6.89) 

2.50  (17.2) 

GLT' 

, 106  psi  (GN/m2) 

0.7  (4.82) 

1.20  (8.26) 

VLT 

0.2 

0.25 

stacking  layers  at  different  angles  of 
inclination  to  the  axis  of  the  joint.  The 
elastic  properties  of  such  a joint  were 
estimated  by  doing  the  usual  property  trans 
formations.8  A simple  laminate  theory  was 
used  to  estimate  the  average  elastic  prop- 
erties of  the  laminate.  This  simple  theory 
ignores  the  out-of-plane  phenomena  such  as 
interlaminar  shear  stresses,  etc. 


Figure  2 shows  the  general  distribution  of  stresses  on  the  fastener  hole 
periphery.  In  some  cases  the  distribution  of  ar  has  the  appearance  of  "half  cosine 
distribution.  As  is  true  in  the  anisotropic  solutions,  sizable  variations  can 
take  place  in  any  direction  depending  upon  the  configuration  of  material  and  geom- 
etries. An  example  of  interest  is  given  in  Figure  2 for  the  curve  marked  with 
squares,  where  e = 1.5  and  the  material  is  0°  glass  epoxy.  In  Figure  2 where 
e/R  = 1.5,  the  method  of  solution  places  certain  limitations  on  geometry  for  its 
successful  results.  One  limitation  is  due  to  the  convergence  of  a Laurent  series 
for  any  awkward  shape.  For  example,  when  either  e/r  or  s/r  is  close  to  unity,  the 
boundary  conditions  are  satisfied  in  a poor  way.  One  more  observation  suggests 
that  for  large  values  of  e and  s the  solution  gradually  merges  into  the  solution 
for  infinite  plate.  The  method  of  solution  was  relatively  faster  in  convergence 
for  periodic  arrays  than  for  the  single  pin  case  with  the  same  number  of  terms 
needed  for  a given  level  of  accuracy  in  the  satisfaction  of  the  boundary  con- 
ditions. The  bulk  of  the  results  tried  were  for  the  perfect  fit  case,  i.e.,  pin 
radius  equals  the  hole  radius.  This  method  was  demonstrated  to  be  suitable  for 


Figure  2.  Radial  stress  distribution  around  fastener  hole  (D  = 2). 


8.  HALPIN,  J.,  ASHTON,  J.,  and  PETIT,  P.  Primer  on  Composite  Materials  Analysis.  Technomic  Publishing  Co.,  1969. 
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clearance  fit  cases,  particularly  when  the  difference  between  pin  and  hole 
radius  was  small.  One  such  example  is  shown  in  Figure  2 labeled  with  triangles, 
applying  to  a pin  radius  of  0.99974. 

It  is  customary  to  compute  stress  concentration  factors,  which  in  turn  help 
in  estimating  the  stress  peaks  in  comparison  with  the  mean  value  of  stresses. 

Now  we  shall  define  the  following  parameters  which  are  of  practical  interest  to 
the  designers  of  composite  materials. 


Mean  Net  Section  Tension,  a^g  = P/2(s-D/2)t 
Mean  Bearing  Stress,  aD  = P/2Dt 

D 

Mean  Shearout  Stress,  x = P/2et 

so 

Further,  these  mean  stresses  can  be  divided  by  the  applied  load  a 
kNS  = 5NS^aA  = S/(S_D/2) 


A* 


kB  = °B/aA  = 2s/D 


The  following  ratios  defined  as  stress  concentration  factors  are  the  ratio 
of  the  peak  stresses  to  the  corresponding  mean  values. 


KNS  ‘Ve  = 90,  r = R/aA 


Kr  = a ) /a 
B B max  A 


K' 

so 


x ) 

xy  max  A 


The  quantity  °g)max  is  defined  as  the  compressive  axial  stress,  i.e.,  maximum 
compressive  value  of  on  along  the  circular  hole  (r  = R)  in  the  contact  region. 
Similarly,  xxy)max  is  the  maximum  shear  at  r = R.  It  was  observed  that  the  max- 
imum xXy  occurs  between  polar  angles  of  about  40°  to  70°  in  the  region  of  contact. 

Finally,  the  ratios  of  °e) g - 90° , r = R aB)max,  and  xxy)max  to  aA  are  of 
interest.  These  are  denoted  as 


kns  °e^6  = 90,  r = R/aNS 
kB  aB^max/^B 
kso  Txy^max/,Tso‘ 

The  following  relations  are  apparent: 
= 
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K'  = k K . 

SO  SO  SO 

Solutions  were  computed  for  various  geometries.  The  important  parameters 
for  observation  are  e/D  and  s/D.  The  material  was  high-modulus  graphite  and  the 
fiber  orientation  was  O2  ± 45.  The  relevant  elastic  properties  are  shown  in  Table 
1.  The  boundary  condition  was  for  the  multiple  pin  case.  The  effect  of  varying 
the  geometry  is  shown  in  Figures  3 and  4.  Figure  3 shows  the  effect  of  s/D  on  the 
stress  concentration  factors  Kg,  and  Kso.  These  curves  stop  at  the  indicated 

dotted  lines  because  s cannot  be  less  than  D/2.  Figure  3 also  indicates  the  lin- 
ear trends  predicted  by  the  infinite  plate  solution  which  shows  the  values  of  s/D 
for  which  the  finite  plate  solution  merges  with  the  infinite  plane  solution. 

Figure  4 shows  the  variation  of  the  stress  ratios  with  normalized  edge  distance 
(e/d)  corresponding  to  s/D  =1.  It  is  interesting  to  note  that  the  curves  fast 
approach  a constant  level  for  e/D  > 2. 


id  h 


10  r 


s - 


4 - 


4 - 


2 - 

l 


Figure  3.  Stress  ratios  versus 
fastener  spacing,  02±45°  graphite 
epoxy  - - e/D  = 4. 


Figure  4.  Stress  ratios  versus  edge  distance, 
02±45°  graphite  epoxy  - - s/D  = 1. 
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Figure  5 shows  some  variations  of  ax  at  various  cross  sections.  We  have 
observed  that  at  a distance  equal  to  twice  the  diameter  away  from  the  hole  the 
stresses  become  nearly  uniform. 


y 

1 L 1 1 ..1 L_ / 

i i 

-8  -7  -6  -5  -4-3  -2  1 

12  3 4 

Figure  5.  Variation  of  axial  stress  profile  along  x-axis, 
02±45°  graphite  epoxy,  e/D  = 2,  s/D  = 0.75,  L/D  = 4 (multipin). 


CONCLUSIONS 

The  results  show  the  complicated  nature  of  the  strain  distributions  in  com- 
posite joints.  The  Laurent  series  approach  brings  out  the  sensitive  points  in 
the  analytic  model  such  as  the  normal  stress  distribution  around  the  hole  as  well 
as  variation  of  the  contact  angle.  The  normal  stress  ar  is  often  shown  to  be  far 
from  the  usual  cosine  distribution.  The  angle  of  contact  is  peculiar  to  the 
elastic  properties  of  the  material  at  the  aspect  ratio  in  consideration.  Theo- 
retically, it  is  also  dependent  upon  the  load  level,  but  for  the  elastic  properties 
of  the  composites  considered  and  the  magnitude  of  the  load  the  joint  would  be 
called  upon  to  withstand,  the  contact  angle  can  reasonably  be  considered  constant. 
The  method  gives  good  solutions  in  the  range  of  geometries  encountered  in  practical 
engineering  design  situations.  The  solution  scheme  breaks  down  when  the  geometry 
is  too  unsymmetrical  about  the  center  of  the  hole.  When  the  geometry  enlarges 
about  the  center  fairly  uniformly  in  all  directions,  the  solution  merges  with  the 
solution  for  infinite  plate. 
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